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ABSTRACT 

A unique feature of deflagration-to-detonation (DDT) white dwarf explosion models of SNe of type 
la is the presence of a strong shock wave propagating through the outer envelope. We consider the 
early emission expected in such models, which is produced by the expanding shock-heated outer part 
of the ejecta and precedes the emission driven by radioactive decay. We expand on earlier analyses by 
considering the modification of the pre-detonation density profile by the weak-shocks generated during 
the defiagration phase, the time evolution of the opacity, and the deviation of the post-shock equation 
of state from that obtained for radiation pressure domination. A simple analytic model is presented 
and shown to provide an acceptable approximation to the results of ID numerical DDT simulations. 
Our analysis predicts a ^ 10^ s long UV/optical fiash with a luminosity of ^ 1 to ~ 3 x 10'^^ erg s^^. 
Lower luminosity corresponds to faster (turbulent) deflagration velocity. The predicted luminosity 
of the UV fiash is an order of magnitude lower than that of earlier estimates, and is expected to be 
strongly suppressed at i > tdrop ~ 1 hr due to the deviation from pure radiation domination. 
Subject headings: supernovae: general — supernovae: individual (la) — white dwarfs — stars: evolu- 
tion — shock waves 



1. INTRODUCTION 

Type la supernova (SN la) explosions are commonly 
accepted to be driven by unstable thermonuclear 
burning of C/O white dwarfs. A key open question 
is whether the nuclear burning front propagates as a 
defiagration o r as a detonation wave. Both deflagra- 
tion models (llvanova et all 119741 : iNomoto et~all I1976L 
119841: iWooslev k, Weaver] 1 1986f rand delayed-detonation 
models, in which a deflagration wav e spontaneous! 



moaeis, m wmcn a aenagration wav e spontaneously 
transforms into a detonation wave (IWooslevI [l99m 
Khokhlovl [19911: lYamaoka et all ri992t iKh okhlov e t al. 
19931: lArnett fc Livnd Il994allbl: [W oosley & Weaver 
I99I; 'H oflich et all 119951 : iHoflkhfc Khokhlovi il996l : 



Iwamoto et al. I Il999n . are claimed to best fit the 
observational data. 

One consequence of a deflagration to detonation tran- 
sition (DDT) is the presence of a shock wave that erupts 
through the progenitor's surface. The detonation wave 
transforms into an ordinary shock wave, behind which 
thermonuclear burning does not release much energy, 
near the edge o f the star, at densities ~ 10®g/cm^ (e.g. 
iPiro et al]|2010[ ) , where the thickness of the burning layer 
becomes comparable to the scale height. Unique signa- 
tures of the presence of a strong shock are an X-ray out- 
burst, which is expected to be produced when the shock 
breaks out from the stell ar edge (Colgate 1974; Falk 1978; 
iKlein k, Chevalieri[T978l ) . and early UV-optical emission, 
which is produced by the expandin g and cool ing shock 
heated outer shells of the progenitor (jFalkll 19781) and pre- 
cedes the emission powered by radioactive decay. 

The observed properties of the shock breakout 
were investigated using analytic order of magnitude 
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estimates (e.g. iMatzner fc McKei Il999t iKatz et al.l 
I20T0I: IPiro et all I2OIOI: INakar fc Saril 120101) and nu- 

merical calculations for particular progenit ors (e.g. 
pnsman fc Burrows 1992: Blinnikov ct al. 200 0l: lUtrobinI 
!2007; Hoflich fc Schacfer 2009; Tominaga et al.' '2009|; 
/Tolstov 2010; Dessart et al. 2011; Kasen et al. 2011). 
An exact description of the time dependent radiation 
emission for non-relativistic shock breakout from a gen- 
eral progenitor (without an optically thick wind) has 
been recently provided in (jSapir et al.l 1201 It IKatz et aD 
l2011t ). This description is accurate during the "planar 
phase" of the expansion, that is as long as the dis- 
tance traveled by the expanding shells is small com- 
pared to the progenitor's radius. Approximate analytic 
descriptions of the emission during the later "spherical 
phase", when the shells expand to radii much larger 
than the original radiu s of the progenitor, were given 
by se veral authors (e.g. i Chevalier 1992; Waxman et al. 
20071 [Chevalier fc Franssoh l2008t iRabinak fc Waxman 
20ld iNakar fc Sari .2010.) . The recent analysis of 
Rabinak fc WaxmanI ( 2010t ) provides an approximate an- 



alytic description including the effects of opacity varia- 
tions due to recombination, which are of particular im- 
portance for compact progenitors with hydrogen poor en- 
velopes, and a method for inferring the relative extinction 
directly from the UV/0 light curves. 

The early emission from a WD undergoing a DDT ex - 
plosion was recently discussed by INakar fc Saril ()2010( ) , 
who gave g enera l order of magnitude estimates, and by 
IPiro et al.l ()2010[ ). who used a more realistic description 
of the propagation of the detonation wave. In this paper 
we improve on the existing analyses by taking into con- 
sideration three effects that have not been addressed ear- 
lier: the modiflcation of the pre-detonation density pro- 
flle by the weak-shocks generated during the deflagration 
phase, the time evolution of the opacity due to recom- 
bination, and the deviation of the post-shock equation 
of state (EOS) from that obtained assuming the pres- 
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sure to be purely contributed by radiation (neglecting 
the plasma contribution). 

A note is in place here regarding the third effect. Dur- 
ing the "spherical phase" the radius r and density p of 
the radiation emitting shells at some given time t is es- 
sentially determined by the expansion velocity v and by 
the opacity n (r ^ vt, t i^pf), and largely independent 
of the initial radius of the progenitor . The initial post 
shock pressure of the radiation emitting shells is ~ pov^, 
where po ^ (rlR^)^p is the initial density. This implies 
that the initial ratio of radiation pressure, aT'^ ~ Po^^j 
to plasma pressure, ^ PoT / p, is (x p^ oc i?* . Thus, 
while neglecting the plasma pressure is an excellent ap- 
proximation for large (core-collapse SN) progenitors, it 
is not a good approximation for SN la progenitors. As 
we show below, taking into account the deviation from 
pure radiation domination leads to strong suppression of 
the flux at t > idrop ~ 1 hr. 

The properties of the emission during both the planar 
phase and the spherical p hase are not very sensitive to 
the i nitial density profile jlSapir et al.l [20111 : iKatz et all 
120111: IRabinak fc WaxmanI 12010^. We therefore discuss 
first, in § 12.11 the opacity and EOS effects using a simple 
analytic model for the spherical phase emission, assum- 
ing a power-law dependence of the progenitor's density 
on distance from the stellar surface and a simple self- 
similar description of the shock propagation (we include 
in § 12.1.21 also a derivation of tdrop using a more realis- 
tic description). The planar phase, which is not much 
affected by the opacity and EOS effects and lasts for 
< 0.1 s for the compact SN la progenitors, is briefly dis- 
cussed in 12.21 We then discuss, in §[3] and §|31 the effect 
of deviations from the simple description of § 12. 1[ using 
detailed hydrodynamical models of DDT explosions. The 
numerical models are described in § [31 and the emission 
properties are described in § 21 Our results are summa- 
rized and discussed in § [51 

2. SIMPLE ANALYTIC ESTIMATES 

2.1. The spherical phase 

For the analysis of this section we assume 
that the shock accelerates down the density 

gradient following the self-si milar solut ion of 

iGandel'Man fc Frank-Kamenetskiil ()1956[ ): iSakurail 
(|1960D . V oc with B = 0.19. This description of the 
shock velocity is valid for explosions i n which the energy 
is rel eased at small radii, r <^ (see lMatzner fc McKed 
Il999f) . i.e. ignoring the fact that part of the energy 
is released by nuclear burning at large radii. The 
distributed energy release is accurately described in our 
numerical analysis (§ [31 and §15). As we show there, 
the simple model described in this section provides 
an acceptable description of the velocity profile of the 
ejecta, and a rough description of the pressure profile 
in regions where radiation dominates the pressure (see 
fig.p. 

The early, t < I day, emission from SNe is dominated 
by the outer < 10"'^ M© shells of the progenitor, which 
are heated by the SN shock and are emitting radiation 
while expanding and cooling. We discuss the opacity 
and EOS effects in § 12.1. Hand ^ 12. 1.2 | respectively , using 
the analytic model of Rabinak fc WaxmanI ()2010D . This 
simple model's assumptions are that the initial density 



profile of the progenitor is given by a power-law of the 
form pq oc S" where 6 = {R^/ri — 1) and is the initial 
radius of the shell, that the post-shock energy density is 
dominated by radiation, and that the post-shock expan- 
sion is adiabatic. Since the results are not sensitive to 
the value of n, we adopt below n — i (the value obtained 
for non degenerate radiative envelopes). 

The model assumptions are similar to those 
used in oth er analyt ic stud ie s of the early emis- 
sion (e.g. 'Chevalier 'Waxman et ah' '2007; 
Chevalier fc Fransson 2008: Rabinak fc Waxman 2010; 



Naka r fc Sanll2010f ). We use the iRabinak fc WaxmanI 
( 20 id) analysis since it includes a realistic descrip- 
tion of the opacity (beyond the Thomson opacity of 
fully ionized plasma used in other analyses). An- 
other comment is in place here regarding the various 
analyses. For the calculation of the luminosity, the 
dif fusion of radiation be l ow th e photosphere is ignored 
in IRabinak fc WaxmanI (l2010fl while it is taken into 
account in iChevalier fc Fran sson' (' 200l: iNaka r fc Saril 
()2010[) . As pointed out in IRabinak fc Waxma n (20101. 
the effects of the diffusion on the luminosity are in- 
deed negligible, and the results obtained with and 

without diffusion are nea r ly identical (e .g. compare 

IRabinak fc WaxmanI IMot INakar fc Sari l20lol) . The 
effects of diffusion are, on the other hand, important for 
determining the spectrum, or color temperature. The 
color temperatures obtained in [Rab inak fc WaxmanI 
(pOTol ) differ from those obtained iu mkar fc Saril ()2010f ) 
due to the more accurate description of the opacity. 

2.1.1. Opacity 

As long as photospheric temperature is well above the 
recombination temperature, the ejecta is nearly fully ion- 
ized and the opacity is time independent and dominated 
by electron scattering. At this stage the bolometric lu- 
minosity and the effe ctive temperature are given by (see 
Eqs. (13) and (15) in IRabinak fc WaxmanI [20Tol) 



L{t) = 3.2 X 10 



and 



39 



pO.85 p 



/0-16(M/1.4MQ)0-69K0;f ^ 



ergs-\ 



Teff(i) =3.5f 



1/4 
'8.5 



0.022 



eV. 



(1) 



(2) 



Here t = 10^i2 s is the time after breakout, i?* = 
10^-^i?8.5 cm, n — 0.2^0.2 cm^ g~^, M is the ejecta mass, 
E = 10^^i?5i erg is the explosion energy and fp is a di- 
mensionless factor that depends on the envelope density 
profile (see Rabinak fc Waxman 2010). K0.2 ~ 1 -t- AT, 
where X is the H mass fraction. The color temperature 
Tcoi is set by photons that can diffuse and reach the pho- 
tosphere after thermalizing in an inner layer. As a result, 
Tcoi is higher than the effective te mperature and is ap- 
proxi mately given by Tcoi ~ 1.2Teff (|Rabinak fc WaxmanI 
[201^. 

When the photospheric temperature approaches the re- 
combination temperature of the ejecta, the approxima- 
tion of constant opacity is no longer valid and the decline 
of the opacity, which is composition dependent, has to be 
taken into consideration. For H poor envelopes, the de- 
cline in the opacity becomes important for T^s ^ 3 eV. 
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The bolometric luminosity for He and C/0 envelopes 
is given in this regime b y (see Eqs. (25) and (29) in 
iRabinak fc Waxman|[20Tol ) 

L[-1(0 = 3.1 X 10^ %o,,5£?;!\n.. ^3°" ergs- 



/o-i5(M/1.4Mo)0-67 3 



and 



(3) 



L[c^/oi(i) =4.4 X 1039 — 



ergs-\ 



/0-i4(M/1.4Afo)°-^^ ^ 

(4) 

respectively (note, that these equation correct for 
a typo that appear ed i n Eqs. (25) and (29) in 
IRabinak fc Wa xman (2 01CII ). where the subscript of R 
should be "12" and not "13"). Here, t = IOH3 s. The 
ef fective temperature is given by (see Eqs. (23) and (27) 
in IRabinak fc Waxmanll2010[) 



Tif'(t) = 1.5 eV/-"-°2i?°Jt3-°-3«, (5) 



and 



-0.017 p0.19j.-0. 35 

/To c to 



(6) 



both with very weak dependence on E and M. 

The composition of the WD's outer mass shells 
at the onset of the explosion is uncertain. The 
estimated m aximum mass of H a nd He shells are 
^ IQ-'^Mrr, (IShen fc Bildsten' '2009a') a nd < IQ-^Mp, 
\lhen fc Tutukovl 11989: Shcn & Bildsten 2009H ) respec- 
tively. As the diffusion sphere propagates inwards to 
larger mass shells it may pass through shells with differ- 
ent compositions. For ^2 ^ 1, the plasma in the photo- 
sphere is fully ionized for all types of compositions, and 
the luminosity and temperature are given by eqs. ([l} and 
^ respectively. For ^2 > 3 the photosphere propagates 
beyond the outer ^ 1O~^M0 and the effective tempera- 
ture drops below 3 eV. At ^2 > 3 we therefore expect the 
radiation to be emitted from shells dominated by He or 
C/0, with luminosity and temperature given by eqs. ([3]) 
and ([S]) or ^ and ([5]), depending on the composition. 
For 1 < ^2 ^3, the detailed properties of the emission 
depend on the amount of H {X as function of mass). A 
comparison of eqs. ([T]) and ([3]) or (|4]) implies that the 
modification of the opacity with time leads to a nearly 
time independent luminosity at t larger than few hundred 
seconds. 

2.1.2. EOS 

As long as the plasma pressure, p, is dominated by 
radiation, and the evolution is adiabatic, the radiation 
pressure of a given fluid element drops like p = Piad cc 



oc r 



pi/3 

density is dominated by the plasma thermal energy we 



For shells at which the thermal energy 



have p = pgas oc p' 



,5/3 



which implies, for radia- 



tion in thermal equilibrium, p^ad p^/^ g<; ^.-8^ js^^ ^-j-^g 
late stages in which we are interested the radiation may 
not be in thermal equilibrium, but is expected to be in 
Compton equilibrium (the number of collisions required 
to modify the photon energy is ~ c/v while the num- 
ber of collisions a photon undergoes over a dynamical 
time is ~ tc/v). For Compton equilibrium, the radia- 



,5/3 



OC r 



tion energy density drops like Prad oc Pgas oc p' 
Thus, at times at which radiation is emitted from shells 
at which the thermal energy density is no longer domi- 
nated by radiation, the radiation energy density in the 



emitting shells and the luminosity are expected to be 
suppressed by a factor ^ (p/po)^^'^ ^ (r/i?,)^-^ to 
~ {p/ PqY^^ ~ (r/i?*)~^ compared to the luminosity ex- 
pected assuming radiation domination (i.e. compared to 
the luminosity given in § I2.1.ip . In what follows we es- 
timate the time <drop, defined as the time at which the 
diffusion front reaches layers in which the ideal gas pres- 
sure is equal to the radiation pressure. For t > idrop a 
strong suppression of the luminosity is expected. 

Assuming that the post shock pressure is dominated by 
radiation (prad — where T is the temperature and 

a is the radiation constant) and that the plasma pressure 
is given by Pgas = PshT/fiirip (where nip is the proton 
mass and p, is the molecular weight), and using the self- 
similar shock de scription (assuming an adi abatic index 
of 7 = 4/3, e.g. IRabinak fc Waxmanll2010l , for details), 
the post-shock ratio of radiation to gas pressure is given 
by 



Prad 



1.2 X 10^Eli'^{p/2) 



Pgas {PO/ g Cm-3)0-53i?0;f4(M/1.4Mo)0-' 



(7) 

(electron degeneracy pressure and pair-production, 
which where neglected in this analysis, do not change this 
result substantially). The post-shock pressure is thus not 
dominated by radiation for 



po > 6.6 X 10' 



1.9 



rI:I{m/iaMq)o-s9 



g cm 



(8) 



Let us next estimate the evolution of the radiation to 
gas pressure ratio, 77, with time. Assuming adiabatic 
evolution of an ideal gas in equilibrium with radiation, 
the density p for which 77 = 1 for a shell with post shock 
density psh = [(7+l)/(7^1)]po and initial pressure ratio 
77 = 77sh, is 

P/P^h = r7,-ie-«("--i). (9) 

To estimate the time when the diffusion sphere reaches 
a layer with i] ~ 1, we use the (time dependent) 
density and pressur e prof iles of the ejecta given in 
IRabinak fc WaxmanI (|2010D . and use the method de- 
scribed there for determining the (time dependent) mass 
of the shell reached by the diffusion front. For time in- 
dependent opacity we find 



rop 



5.3 X 10 



/O-"(Af/1.4M0)O-69 



s, 



(10) 



while for C/0 opacity we find (using the power- 
law approximation for the opacity suggested in 
IRabinak fc WaxmanI [2QToI )) 



idrop « 3 X 10^ 



1.2 



Elf^Rs.5W2) 
/O-i5(Af/l.4M0)O-56' 



(11) 



(in both cases we neglected logarithmic corrections). 

It is illustrative to compare idrop obtained above to 
that obtained using the m ore det ailed description of the 
ejecta profiles suggested bv lPiro et al. (2010). These pro- 
files are obtained under the following assumptions: the 
pre-detonation envelope is in hydrostatic equilibrium and 
the pre ssure is dominated by degeneracy pressure (se e 
eq. IA2p , the velocity of the shock is given by eq. (|A5p , 
and the terminal velocity of a shell is vj = 2vs (instead 
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of f / = (6/7)fs taken by iPiro et al.ll201d . which IS more 
appropriate for the planar phase), the diffusion front is 
located at an optical depth r = c/w/. The pressure in the 
ejecta is calculated using the relation p = psh {p/psh)^ , 
where is the initial pos t shock pressure (correcting 
Eq. (17) of lPiro et al.][2Eol which uses po instead of psh 
in the denominator). For these profiles we find 



Pgas {Po/ g cm 3)0.52 ' 



(12) 



where = lO^wg cm sec~^, and prun — lO^pe g cm^'^ 
are the velocity and density where the detonation wave 
transforms to a shock wave respectively. Thus, the as- 
sumption of radiation pressure dominance breaks for 



po>5xl0\rpr'{p/2y-'gcm-\ 
For constant opacity we find 



(13) 



id: 



Rhvl''4ipl''{p/2)^' 

rop ~ o.u ^ xu (M/1.4Mo)0-5 



3.6 X 10^- 



s, (14) 



while for C/0 opacity we find 



rop 



nl.74 1.76„0.31/,,/o\1.5 

(M/1.4M0)"-44 
2.2. Breakout flash 



(15) 



When the shock reaches layers with optical depth 
T = SMbo'^/{4:^R*) comparable to c/vs, where SMqq is 
the mass exterior to the shell, photons outrun the shock 
and escape, producing a breakout flash. Our simulations' 
pre-detonation profi les (described in jj lA.ip are different 
from those used bv IPiro et al.l (j201C0 . However, as ex- 
pected, the breakout energy is not sensitive to the details 
of the profiles. Using the simulation based (extrapolated) 
profiles derived in §[Bl and the shock description given in 
ii lA.2[ the shock velocity at t = c/vs is mildly rclativ istic 
(for Ko.2 = 1)- For mild-relativistic shocks, Vg in Eq. (jASP 
is replaced with TsVs, where T g = [1 — {vs/c)'^]~^' 
(|Tan et al.l 120011: IPiro et al.ll2010D . and we find that the 
breakout energy is 



BO 



10^ 



n2.a„,1.16 fU.05/'9^ /,,V>:21 

(A//1.4Mo)0-i6(r,/2.1)i-i6K<]-84 



erg. (16) 



i?BO depends weakly on fp = — /3), where 

1 — 13 — L/Ledd, the ratio of the luminosity escaping 
the pre-detonation progenitor to the Eddington lumi- 
nosity, is not accurately determined by our simulations. 
Our results is in agreement with prev i ous estimates 
(llmshennik et al.l 119811 : IPiro et al.l l2010t iNakar fc Sari 
|201CI[ ), su ggesting that the h i gher breakout energy ob- 
tained bv lHoflich fc Schaefeii ()2009[ ) is due to the coarse 
numerical resolution near the stellar edge. The breakout 
fiash is spread over i?»/c (the photon diffusion time at 
shock breakout, which is smaller and is not weakly de- 
pendent on fp, has little effect on the observed proper - 
ties of the emissio n) . As shown bv lBudnik et al.l ()2010D : 
iKatz et all ()2010D , for mildly relativistic breakout veloc- 
ity (v/c > 0.2), a non-thermal spectrum extending to few 
hundred keV is expected. 



3. NUMERICAL CALCULATIONS 

We carried out one dimensional (ID) simulations of 
DDT SN la explosions, in order to obtain the outer 
ejecta profiles (i.e. density, pressure and temperature 
profiles) , that determine the emission of radiation during 
the "spherical phase". Some aspects of burning, such as 
fiame instabilities and turbulence, cannot be addressed 
in ID simulations and are phenomenologically treated by 
using a parameterized defiagration velocity. Despite this 
caveat, we expect the simulations to describe the outer 
ejecta profiles truthfully. 

The simulations hav e been perfo rmed using the Vul- 
can/lD (VID) code (|Livnel Il993f ). that incorporates 
Lagrangian hydrodynamics, general EOS and nuclear 
burning. The EOS includes contributions from fully 
ionized gas and radiation at LTE. Free electron pres- 
sure and energy are computed using tabulated val- 
ues of partially degenerate electron-positron gas. For 
burning we used the Alpha network of 13 elements 
from "'He to ^^Ni. The nuclear reaction rates are 
ta ken from the NON-SMOKER database as described 
bv lRauscher fc ThielemannI (|2001[ ). 

We assumed the " standard" scenario of thermonuclear 
runaway in a Carbon-Oxygen WD that approaches the 
Chandrasekhar mass by accretion from a companion. We 
assumed equal mass fractions for the two species, ^^C 
and ^^O. For a given central density and temperature we 
integrated the equation of hydrostatic equilibrium out- 
wards under the assumption of constant entropy. We 
then modified the central density until the desired mass 
is achieved. Using variable zone masses we obtained a 
reasonable smooth mass distribution which goes down 
to ~ 1O~^M0 in the outer cells. The total number of 
zones used was roughly 1000. The number of zones was 
increased to over 4000 in one simulation to test for con- 
vergence. 

Our simulations start with a defiagration phase where 
the defiagration front propagates from the center of 
the progenitor outwards, with a parameteriz ed local 
speed of the form used by iHoflich et al.l (|2002t ) , where 
the deflagration velocity is take as the maximum be- 
t ween the Laminar de flagration velocity (see Eq. (43) 
of iTimmes fc Wooslevl[l99 2') and turbulent velocity. For 
the turbul ence velocity we use the para metrization sug- 
gested bv iDominguez fc HoflichI (pOOOl ). /turb(5^i)^/^ 
where A is the Atwood number, g is the gravitational 
acceleration, L = {dr/dp)p is the characteristic pressure 
length scale and /turb = 0.2 and 0.4 for "slow" and "fast" 
simulations respectively. The values for ftnrh were cali- 
brated from 3D calculations (Khokhl ov et al.l l2002\ 

During the deflagration phase, the progenitor expands 
and its outer layers move with a typical velocity of 3 to 
6 X 10* cm s~^ depending on the defiagration velocity. In 
regions where the expansion velocity is higher than the 
matter sound speed, weak shock waves form and heat the 
plasma. This leads to modifications of the density and 
pres sure profiles of the progenitor, as discussed in detail 
in § lA.ll The fully ionized gas EOS is not valid at the 
outer parts of the progenitor, leading to inaccurate initial 
profiles in these regions. These profiles, however, are 
modified by the heating during the defiagration phase, 
and therefore have little effect on the final ejecta profiles. 

The deflagration front propagates until it reaches a 
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critical density pcrit, where we induce a detonation wave 
by increasing the deflagration speed to nearly sonic over 
few zones. Three values are taken for pcrit^ 1, 2 and 
2.5x10^ g cm"^. These are denoted by 1, 2 and 2.5 in 
the simulation names. The detonation front propagates 
outwards until the entropy produced by the burning of 
elements becomes smaller than the entropy that is pro- 
duced by the compression, at which point the detonation 
transforms into a shock wave. The transition of the det- 
onation to a shock wave as well as the shock propagation 
are discussed in detail in § IA.2I We follow the expansion 
of the progenitor after the shock wave reaches the stellar 
edge, neglecting the effects of radiation transport. 



TABLE 1 
Simulation Parameters 
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Name" 


Color 


'dct 




7p^ 




Ell 






[s] 






[10*cm] 




DDTls 


blue 


3.05 


3.5 


0.95 


5.3 


1.27 


DDT2s 


green 


2.8 


7.1 


0.89 


4.3 


1.52 


DDT2.5S 


red 


2.75 


7.9 


0.89 


4.1 


1.55 


DDTlf 


cyan 


1.75 


1.3 


1 


4.9 


1.31 


DDT2f 


magenta 


1.6 


6.7 


0.86 


3.9 


1.55 


DDT2.5f 


black 


1.6 


6.7 


0.86 


3.9 


1.58 



"1,2 and 2.5 correspond to DDT transition at Pcrit = 

10^ g cm~^x 1, 2 and 2.5. "s" and "f stand for slow and 

fast deflagration velocities (cf. §[3ll 
^ The time at which detonation is induced. 

Fit parameters for Po{po) a^t ts = idet a-nd 

po < 10* g cm-3 (eq. lIMt ). 



Progenitor radius at t = t^ct ■ 
' Total ejecta energy in units of 10^^ 



erg. 



4. UV/OPTICAL EMISSION DURING THE 
SPHERICAL PHASE 

4.1. Density and pressure profiles of the outer parts of 

the ejecta 

We derive the emission of radiation using the density 
and temperature profiles at the end of the simulations, 
ts — 100 s (note that tg is measured from the onset of 
deflagration, while we use t to denote the time measured 
from breakout). At this stage the ejecta has expanded 
significantly and the radius r of a relevant mass shell 
is ^ i?, and is approximately given by tvf, where the 
shell's velocity w/ is approximately its terminal veloc- 
ity. The simulations' mass resolution at the outer cells, 
~ 3 X 10~^Mq, sets a lower limit to the time from which 
the emission can be reliably calculated based on the sim- 
ulations, t >3 X lO^s. 

We neglect the reheating by the radioactive decay of el- 
ements synthesized in the explosion process (valid for the 
outer layers of the ejecta, in which only a small fraction 
of the elements are fused) , and extrapolate the profiles to 
later times assuming that the ejecta expand adiabatically 
and that the fluid elements reached their final velocity. 
We use an EOS of radiation in equilibrium with an ideal 
gas (/X = 1.7455 for fully ionized equal mass fractions 
of C and O). Note that although we assume full ioniza- 
tion for the EOS, we do not make this assumption for 
the opacity. This approximation for the EOS is accurate 
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Fig. 1. — Pressure and velocity as function of electron-scattering 
optical depth obtained from the simulations (solid lines) at ts = 
100s, compared to the a nalytic es timates of Rabinak & Waxma3 
(2010) (dashed lines) and'Piro et al. (201 0) (d ash-dotted, corrected 
as explained in the text preceding Eq. ||12| |). Line colors corre- 
spond to different simulation parameters as given in Table [T] "x" 
marks indicate the position where the radiation and gas pressure 
are equal. 

as long as the thermal energy density is dominated by 
radiation, and becomes less accurate when the plasma 
pressure becomes significant. This implies that our light 
curves are not accurate at i > tdrop- 

In Fig. [T] we compare the pressure and velocity pro- 
files at th e end of the simulat i ons to the profiles sug- 
gested by iRabinak fc WaxmanI (|2010D (for n = 3) and 
iPiro et ahl (2010 corrected as explained in the text pre- 
ceding Eq. (fT2|) ). denoted by RW and PCW subscr i pts re - 
spectively. For the model of IRabinak fc WaxmanI ()2010l ) 
we used £^51 = l,i?8.5 = l,andM = 1.4Mq, similar 
to the values obtained from the simulations (see table 
[J). The discontinuity in the PCW profiles results from a 
transition between regions with different pre-detonation 
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density profiles. Both the RW and PCW models do not 
properly describe the pressure behavior at large optical 
depth, where the pressure is no longer purely dominated 
by radiation and thus drops faster with time, as explained 
in§[2T2l 

For the lower deflagration velocity the pressure is 
higher by ~ 70% at 10^ < Tes ^ 10^, where Tes is the op- 
tical depth for electron scattering (ko.2 = !)■ This differ- 
ence is due to the fact that for the low velocity explosions 
the envelope is heated to higher temperatures during the 
deflagration phase. Thus, for a given pre-detonation den- 
sity po the mass fraction 6M (integrated mass to the stel- 
lar surface) is larger compared to that in high velocity 
explosions. For 7 = 4/3, as Tes oc SM{tVs)~^ and the 
density in the ejecta p oc 5M{tVs)~^ , the pressure in the 
ejecta p oc v-^{p/po)'^ oc {SM/ poY^'^Tcst^^ , implying that 
for lower po at given 6M the pressure is larger at a given 

4.2. L and T^ff. 

We use the density and pressure profiles obtained as 
described in § 14.11 to derive the properties of the emit- 
ted radiation. Due to the deviation from pure radiation 
pressure domination, the pressure profiles deviate signif- 
icantly from p (X T (which is valid for the core collapse 
progenitors). This implies that photon diffusion may 
have a non negligible effect on the predicted luminosity 
(jRabinak fc WaxmanI 12010} ). We therefore estimate the 
luminosity at time t as dE/dt, where E[t) is the radiation 
energy that escapes by diffusion up to time t. The diffu- 
sion depth is determined as the depth at which r — c/vf, 
and T is calculate d with the opacity taken from the OP 
project tables f|Seatonll20Q5[ ). The effective temperature 
is calculated as Tee = {L/iTrr'^^aY/^, where rph is the 
radius of the photosphere (r = 1). The temperature at 
the photosphere (t = 1) is typically larger than Toff by 
10% 20%, which implies that neglecting diffusion leads 
to an overestimate of L by a factor of 1.5-^2. The simula- 
tions' resolution was not fine enough to allow an accurate 
determination of the color temperature. 

The resulting luminosities and effective temperatures 
are given in Fig. [51 The luminosity depends weakly on 
Pcrit, and is lower for larger deflagration velocities, re- 
flecting the lower pressure at given t^s (see Fig.[T]). Com- 
paring the results presented in the figure with eqs. ([6]) 
and PT|) . we find that the simple model of 12. II provides 
a good description of the properties of the emitted radia- 
tion. For the slow deflagration velocity the luminosity is 
similar to that predicted by the simple analytic analysis, 
while for the larger velocity it is 2-3 times lower. 

An estimate of the luminosity expected for a pure He 
composition of the outer shells may be obtained by ex- 
trapolating the simulation profiles assuming an EOS with 
p = 2, and calculating the emission using the opacity of 
He. The luminosity obtained in this case is ~ 50% lower 
and the effective temperature is ~ 15% lower than the re- 
sults obtained for C/0, consistent with the simple model 
of §11111 

5. SUMMARY AND DISCUSSION 

We have presented in § 12.11 a simple model of the 
UV/optical emission during the spherical phase, assum- 
ing a power-law density profile of the progenitor that is 
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Fig. 2. — Bolometric luminosity and effective temperature de- 
rived from the simulations. Line colors correspond to different 
simulation parameters as given in Table[T] t = islO^s. 



not modified during the defiagration phase, and a self- 
similar shock propagation, ignoring the fact that part 
of the explosion energy is released by nuclear burn- 
ing at large radii. The predicted luminosity and effec- 
tive temperature are given in eqs. (ITl)-([6]) for different 
compositions of the outer shells of the progenitor. For 
t > 300 s the photosphere propagates beyond the outer 
~ lO'^ AffT^, the estimated maxi mum mass of a possible 
H sheh UShen fc Bildstenll2009a[) . and the effective tem- 
perature drops below 3 eV. At this time we therefore ex- 
pect the radiation to be emitted from shells dominated 
by He or C/0, with luminosity and temperature given 
by eqs. ([3]) and jS]) or (|4|) and ([6]) respectively. At ear- 
lier times the detailed properties of the emission depend 
on the mass fraction of H. The simple model predicts a 
nearly time independent luminosity, L ~ 10^^'^erg/s, at 
300 s < t < 10^ s, and a strong suppression of the fiux at 
t > idrop ~ 1 hr (see eas.fT4land[T5l) due to the deviation 
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from pure radiation domination. 

Comparing the results of the model to those of detailed 
ID simulations, we find that it provides an acceptable de- 
scription of the velocity profile of the ejecta, and a rough 
description of the pressure profile in regions where radia- 
tion dominates the pressure (see fig. [TJ a detailed discus- 
sion of the pre-detonation profile, which is significantly 
affected at po ^ 3 x 10'' g cm~'^ by weak shocks during 
the deflagration phase, is given in the appendix). Thus, 
we expect the model to provide an acceptable description 
of the luminosity at t < idrop ■ 

Indeed, we find that eqs. (j?]), © and ([TT|) provide an 
acceptable description of the properties of the emitted 
radiation. We find that the luminosity depends weakly 
on pcrit, the density at which detonation is initiated, and 
is lower for larger deflagration velocities. For the slow 
(turbulent) deflagration velocity (see §[3]for a description 



of the turbulent velocity prescription) the luminosity is 
similar to that predicted by the simple analytic analysis, 
while for the larger velocity it is 2-3 times lower (see 
fig.©- 

As expected, we find that the breakout emission dur- 
ing the planar phase is not sensitive to the density pro- 
file modifications. We find Ebo ~ 10'"' erg (see eg. ITC)). 
spre ad over i?^,/c , in ag r eement wit h previous esti mates 
(Ims hennik et al.l 119811: iPiro et al.l ^ 
2010). As shown bv I Budnik et"al 



2010; Nakar fc Sari 
(|201Q ): iKatz et al 



(2010), for the predicted mildly relativistic breakout ve- 
locity {v/c > 0.2), a non-thermal spectrum extending to 
few hundred keV is expected. 

We thank A. Glasner and D. Kushnir for useful 
discussions. This research was supported in part by ISF, 
UPBC, GIF and Minerva grants. 



APPENDIX 

A. ANALYSIS OF NUMERICAL SIMULATIONS' RESULTS 

A.l. Pre-detonation pressure and density profiles 

We consider the pressure and density profiles of the outer layers of the progenitor, which are modified by weak 
shocks generated during the deflagration phase, focusing on the proflles at t = tdot- During the deflagration phase the 
progenitor WD expands. In regions, where the expansion velocity is similar to or higher than the sound speed, weak 
shock waves form, which heat the envelope and modify its pressure profile. During this phase, two physical processes 
affect the pressure profile: shock wave heating and diffusion of energy. Diffusion affects significantly regions in which 
the photon diffusion time (to the stellar surface) is similar to the deflagration time. The outer part of the progenitor's 
envelope may be divided into three regions distinguished by the physical processes which affect them: (1) The inner 
region, which is negligibly affected by shock heating and diffusion; (2) The intermediate region, which is affected by 
shock heating but is negligibly affected by diffusion; (3) the outer region, which is affected by both shock heating and 
photon diffusion. The mass resolution in our simulations does not allow us to fully resolve region (2). We denote by 
(2A) and (2B) the parts of this region which are resolved and unresolved by the simulations respectively. Below we 
describe the simulations' results for regions regions (1) and (2A). The unresolved regions, (2B) and (3), are discussed 
in §11 

Let us first discuss the pressure profile expected at the outer shells of the progenitor based on simple theoretical 
arguments. For simplicity, we neglect the shells' self-gravity and thickness (relative to the progenitor radius at the 
onset of detonation). We denote the density of a shell at the onset of detonation by po- Assuming that the star is in 
hydrostatic equilibrium, we have for the outer layers 

SM{po) « , (Al) 

where Al is the total mass of the progenitor, Pq is the pressure, and SM is the integrated mass from the point where 
the density equals po out to the stellar surface. For a polytropic EOS, Pq oc pj^^^", the density profile is oc (5", where 
S={R,/r-l). 

We consider first the inner region (1). In this region the pressure is negligibly affected by the deflagration process 
and is therefore dominated by the electron degeneracy pressure. In this case Po(po) is polytropic with index n = 3/2, 
and may be approximated by (following iPiro et al.ii2010 ') 

5/3 

il2 / PO \ , „_-3 



Pdcs = 9.91 X 10'^ ^ ^ erg cm'^ (A2) 

VMe g cm V 

where fie is the molecular weight per electron. The sound speed in these layers is Cs = 4.1 x 
106(po/ gcm-3)i/3^^/^ cm s~\ 

Let us consider next region (2A). The simulations show that when the deflagration front reaches a density pcrit ~ 
10'' g cm^'^, at time idot, the outer layers of the progenitor expand at a speed 3 6 x 10* cm s^^. The pressure profile 
deviates from that of region (1) at densities ~ 3 x 10^ g cm^"^, beyond which it becomes flatter and nearly linear in po. 
This change in the density profile is the result of shock heating during expansion. We use the simulations to obtain 
an approximation for Po(po)- In Table [1] we give the best fit parameters for a pressure profile of the form 

Po(po) = X 10^5 ( ^° X ' erg cm~3, (A3) 
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Fig. A1. — Left: P o(po ) obtained in the simulations at t = t(jet (solid lines) compared to the approximations of Eq. IIA2I I (black dash 
dotted line), and Eq. | |A3|| (with the values given in Table [l] dashed lines). The pressure is normalized to the pressure given by Eq. ||A4| | 
for TiokcV = 1 and fi = 2. Right: 5M(po) obtained in the s imul ations at t = fdct (solid lines) compared to the approximation of Eq. llAlll . 
with Pq given by Eq. I IA2II (black dash dotted line) or Eq. I IA3II (dashed lines). The masses are normalized to the mass given by Eq. IIAll l 
with Pq = Pthrm (with TiokcV = 1 a-iid M = 2). Different colors correspond to different simulation parameters as specified in Table [l] 



for t — tdot and po < 10'' g cm^'^. In all of t he fits the index 7„ is in t he range 0.9 1, implying a nearly isothermal 
profile, in agreement with the simulations of 'Hofl ich fc Schaefed (|2009l see figure 1 in that manuscript), which have a 
similar resolution to that of our simulations, < 1O~^M0. 

The pressure in region (2 A) is different than that assumed bv llmshennik et alj ()1981[ ): |Piro "eFaI1(l20l0l ): INakar fc Safil 
((Moh (for the same po range). In particular, it corresponds to a negative polytropic index {n < 0) and does not 
correspond to a density profile declining as a power of d (distance from the progenitor's edge). The density profile 
we find falls approximately exponentially with r, similar to the profile produced by an isothermal gas. The shells of 
region (2A) determine th e em ission on time scales of minutes to hours. 

In the left panel of Fig. lA.ll we compare the results of the simulations to the approximations of Eqs. (|A2|) and (|A3|1 . 
Since the profile is nearly isothermal at the outer regions, we normalize Pq to 

Pth„n(po) ^ 4.79 X 10^^ (^^) erg cm"^, (A4) 

V g cm -^J p,/2 

where T — lOTiokcV keV and fi is the molecular weight. The approximation of Eq. (|A2|) is good to a few percent in 
the density range 3 x 10'' g cm^'^ < po < 10^ g cm~^, and the approximation of Eg. (|A3[) is good to tens of percent 

for pq < 10^ g cm^'^. In the latter density range, Po is given by Pthrm to within a factor of 2. 

In the right panel of Fig. IA.l|we c ompa re S M{po) obtained in the simulations to th at g iven by Eq. (jAip using the 
pressure approximations of Eqs. (jA2p and (jA3l) (with the values given in Table [1]). Eq. (|A1[) with Pq given by Eq. (|A2p 
is accurate to a factor ~ 2 in the density range 3 x 10^ g cm~^ < po < 10^ g cm^'^. The discrepancy is partially due 
to deviations from hydrostatic equilibrium and part ially due to the fact that R„ over estimates r by up to 15% in 
this density range. Eq. (jAip . with Pq given by Eq. (jA3[) . is accurate to a few tens of a percent in the density range 
Po < 10^ g cm~^. In this density range, substituting Pthrm for Pq gives an approximation for SM, which is good to up 
to a factor of 2. 

A. 2. Shock propagation 

As the detonation wave propagates in the declining density profile of the star, the entropy produced by the burning of 
elements decline, until eventually this entropy is smaller than the entropy that is produced by the shock compression. 
At that point the detonation wave transforms to a shock wave. The shock wave continues to accelerate in the declining 
density profile. In regions where the densi ty falls as i5", shock acceleration is described by the G andel'Man-Frank- 
Kamenetskii-Sakurai self similar solutions (iGandel'Man fc Frank-Kamenetskiilll956l ; ISakurailll96ClD . In regions where 
the de nsity falls exponentially with 5, the acceleration is described by the self similar solution of iGrover fc Hardyl 
(jl96(3[) . In both types of solutions, the shock velocity is given by 

Vs^V,^n(—] (A5) 

VPrmi/ 
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Fig. A2. — Left: Shock velocity (blue) and te rmin al velocity {ts = 100 s, red) obtained in one simulation (with higher dumping resolution), 
normalized by the shock velocity given by En. ||A5l l. with /3i = 0.18. Right: Terminal velocities {ts = 100 s, red) obtained in the simulations, 
normalized to the approximation of Eq. l)A6| l with /3i = 0.18 and /„ = 1. Different colors correspond to different simulation parameters as 
specified in Table [T] 



For a radiation dominated shock, with post shock adiabatic index of j<, = 4 /3, (3i w 0.19 for density profiles dechning 
as a power law with n = 3,3/ 2 (jGrassberd 119811: iMatzner fc McKed 119990 . and /3i w 0.176 for exponential density 
profiles (jGrover fc Hardvlll966[ ). The self-similar solutions also show that the terminal velocity of the shells, vj, is 
oc Vs- For power-law density profiles Vf = 2vs, while for exponential densit y profiles v t = l.S f.s. 

The transition from a detonation wave to a shock wave was analyzed bv lPiro et al.l (|2010l ). The transition density 
Prun was defined as the density at which the "induction length" A = (q/e)pdet, where q,e are the energy and energy 
generation rate of the detonation process respectively, is smaller than the scale over which the progenitor density 
changes, assuming pdctis given by the Chapman- Jouguet detonation velocity at the corresponding density. For a WD 
equally composed of ^^C and ^^O they foun d pmn ~ 2 x 10^ g cm^'^, in agreement with our simulations (although the 
equation for e was take in lPiro et al.l (|2010 ) at the wrong electron screening regime.) I t should be noted here that th e 
values of Vrun, Prun may be affe cted by instabilities which are not considered here (e.g. iDommguez fc KhokhlovllMll) . 

In the left panel of Fig. IA2I we compare the shock velocity and the terminal velocity obtained in a simulation to 
the approximation of Eq. (jASp . This fi gure shows that although the density profile changes from a power of S to 
an exponential, the approximation Eq. (jA5l) holds. The figure also shows that the values estimated analytically for 
Piun and Vrun are accurate to ~ 10%. The ratio Vf/vs is not independent of po, and declines from a ratio of w 2.4 
at Po = 10^ g cm~^ to « 1.4 at po = 10 g cm~^. The decline in the ratio Vf/vg with pn can be attributed to th e 
transition from a power-l aw t o an exponential density profile and to spherical affects (see IMatzner fc McKe3ll999| ). 
In the right panel of Fig. IA2I we compare the velocity at tg = 100 s (terminal velocity) to the approximation 

V/ = 2.5/„z;,(po/Prun)"-°'. (A6) 



B. EXTRAPOLATING THE PRE-DETONATION PROFILES TO 5Mbo 

Our simulations do not resolve regions with densities pQ < 10^ g cm~'^ at t^ct and we only consider them qualitatively. 
These regions determine the emission at t < 3 x lO^s (cf. § [2|) and have little effect on the emission at later times. 
Assuming that the pressure profile given by Eq. (|A3|) continues to densities ^ 10^ g cm~^, the pressure becomes 
dominated by radiation at po ~ 1 g cm~^, the flux reaches ^ gc/K and the luminosity approaches the Eddington 
luminosity, L^dd- In th is case, the lumino s ity is L w L^dd — 4:TtcGM/k ~ 10'^^ erg s~^, which is similar to the 
luminosity obtained by iHoflich fc Schaefe? (200i l) befo r e the transition to detonation, and is larger by 3 orders of 
magnitude than the luminosity assumed by Piro et al.l ()2010[ ). The simulations indicate therefore that the heating 
process is efficient in the sense that it leads to L ~ L^dd during the deflagration phase, that this upper bound is 
reached. 

It is reasonable to assumed that near the stellar surface, in regions where the diffusion time to the stellar surface is 
<^ idct, the luminosity is constant. Assuming that the envelope is at hydrostatic equilibrium, that opacity is constant 
and that the pressure is a sum of a radiation pressure and an ideal gas pressure with mean molecular weight p,, we 
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have in this region 

where 1 — P — L/Ledd- The large L obt ained in the simu lations indicates that fi3 = /3~^(1 — /?) ^ 1- This pressure 
profile is similar to the one suggested by iPiro et all ()2010( ) for the outermost layers but its normalization is higher by 
an order of magnitude. For this profile, the diffusion time to stellar edge is 

j-2/3n4 , .5/3 



^ (M/1.4M0)2(/x/2)8/3 (, g cm-3 ) ^^^^ 

where k ~ 0.2ko 2 cm^ g~^ is the opacity. The density for which tdiff' ~ ^dct is roughly po 1 g cm~^ (for tdct given 
in Tabled]). 
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